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Fast quantum control and light-matter interactions at the 10,000 quanta level 
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Fast control of quantum systems is essential 
in order to make use of quantum properties be¬ 
fore they are degraded by decoherence. This 
is important for quantum-enhanced information 
processing [lj, as well as for pushing quantum 
systems into macroscopic regimes at the bound¬ 
ary between quantum and classical physics mm- 
“Bang-bang” control attains the ultimate speed 
limit by making large changes to control fields 
on timescales much faster than the system can 
respond BUS), however these methods are often 
challenging to implement experimentally. Here 
we demonstrate bang-bang control of a trapped- 
ion oscillator using nano-second switching of the 
trapping potentials. We perform controlled dis¬ 
placements which allow us to realize quantum 
states with up to 10,000 quanta of energy. We 
use these displaced states to verify the form of the 
ion-light interaction at high excitations which are 
far outside the usual regime of operation. These 
methods provide new possibilities for quantum- 
state manipulation and generation © along¬ 
side the potential for a significant increase in op¬ 
erational clock speed for ion-trap quantum infor¬ 
mation processing [8l ,9]. 

Quantum control involves the preparation, manipu¬ 
lation and exploitation of high-purity quantum states, 
which constitute a fundamental resource for quantum in¬ 
formation processing [lj. For quantum harmonic oscilla¬ 
tors, experiments have demonstrated such control using 
confined optical and microwave fields mm, and using 
the mechanical oscillations of single trapped ions CMS]. 
Most of this work was performed in the adiabatic and res¬ 
onant regimes, where control timescales are much slower 
than the system’s natural frequencies. “Bang-bang” con¬ 
trol is the opposite extreme, in which timescales are so 
short that changes to the system can be considered to be 
instantaneous mmm- For trapped ions, bang-bang 
control has been performed using ultra-fast laser pulses 
HMD- However, optical fields interact with the electric 
dipole or quadrupole moment of the atom, and the oscil¬ 
lator states prepared in this way have thus far been lim¬ 
ited to excitations of less than 20 quanta. Electric fields 
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FIG. 1: Experimental sequences: a) The ion is first cooled 
close to the ground state of motion in a harmonic trapping 
potential centered at x = 0, and the spin state is pumped to 
|^) (see Supplementary Information). At t = 0 the poten¬ 
tial center is suddenly displaced to xa and the ion’s motional 
state becomes | — <ao). b) Evolution of the motional wave func¬ 
tion in the position-momentum phase-space. After a time At 
in the displaced potential, the trap center is switched back 
to x = 0, leaving the ion in a coherent state of size given 
by equation ([Tj) . c) Sequence of laser-pulses and potential 
displacements for monitoring the time evolution of coherent 
states, d) Sequence for studying the light-matter interaction 
with highly-excited states. In this sequence, At is an integer 
number of times j of the oscillator period. 


acting directly on the charge (electric monopole) can be 
used for achieving a strong-interaction regime. This is 
a key element for transporting trapped-ion qubits in a 
scalable quantum-computing architecture 0 [20], where 
operations have been performed on timescales compara¬ 
ble to the oscillator frequency In these examples 

the speed was limited by the bandwidth of the filtered 
control voltages applied to the electrodes. 

In this Letter, we demonstrate bang-bang control over 
trapped-ion harmonic oscillators, which we use to gener¬ 
ate and characterize quantum states with up to 10,000 
quanta. To achieve this, we wire the trap electrodes ac¬ 
cording to a new scheme: we prepare constant, filtered 
voltages at the inputs of switches placed in vacuum close 
to the trap electrodes, and select the input to be wired 
to the electrode by means of digital control pulses (for 
further details about the experimental setup, see Supple¬ 
mentary Information). For this setup we have measured 
switching times of nano-seconds, which we think were 
limited by our measurement apparatus [8]. This is much 
faster than the ion oscillation period, meaning that we 
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can induce quasi-instantaneous changes to the oscillator 
states (Fig. ©)• Using these techniques, we map out 
the ion-light coupling strength for high excitations, and 
track the oscillation trajectory of a single ion throughout 
multiple oscillation cycles. 

In order to examine the timing quality and repro¬ 
ducibility of our control, we prepare a scenario which 
is highly familiar in classical physics. Starting with a 
particle at rest at the bottom of a harmonic potential, 
we suddenly displace the latter by a distance x&. For an 
instant, the particle finds itself at rest but in an excited 
state, and subsequently oscillates in the new potential. 
In our experiments, this situation is reproduced in the 
quantum regime. In the ideal case, the initial state is 
the quantum ground state |0), which has a minimum- 
uncertainty wave-packet of root-mean-squared extension 
x$. The action of the sudden displacement produces a 
coherent state |—ao) with qlq = x&/(2xo), which oscil¬ 
lates at a frequency cj m while maintaining the form of 
the wave-packet. Switching the potential minimum back 
to x = 0 after a time At leaves the oscillator in a coherent 
state of size 

M = |ck 0 | \/2 [1 - cos(w m A t)\ (1) 

in the original potential (Fig. [l]b). 

The trapped-ion oscillator we consider also has internal 
degrees of freedom, from which we can isolate a two-state 
pseudo-spin system (||), |t)) with transition frequency 
cjo- This transition can be parametrically coupled to the 
oscillator using electromagnetic radiation at frequency 
i^l = + 5u; m . For integer s, resonant transitions can 

be induced on the s th sideband between the states || 
)| n) |t)| n + s ) with a Rabi frequency [23] 

tt n ,n+s = fio \(n + s|e lfea!:E |n)| , (2) 

where Uq is a constant which contains the internal-state 
coupling strength and the electric field amplitude, and 
k x is the size of the radiation wavevector projected on 
the oscillation direction. For k x x o <C 1 and at low exci¬ 
tations, a small parameter expansion of the exponential 
is sufficient to describe the resulting dynamics. This is 
the regime in which most quantum control experiments 
with trapped ions have been operated in the past, and 
as a result coherent control has only been demonstrated 
for | s | < 2. In the work presented below, we map out the 
excitations for 0 < s < 5, confirming the dependence of 
the light-atom interaction strength for excitations up to 
n « 10,000. 

In our experiments, the oscillator is the axial mo¬ 
tion of a single 40 Ca + ion confined in a radio-frequency 
trap with a frequency between u) m / (27r) ~ 2.35 MHz 
and 2.53 MHz, and corresponding values of xq ~ 7 nm. 
We use the dipole-forbidden transition at wavelength 
A « 729 nm to define a pseudo-spin system between 
levels ||) = \L = 0, J = 1/2, Mj = —1/2) and |t 
) m \L = 2, J = 5/2, Mj = —5/2). The wavevector 
of the laser addressing this transition makes an angle 


6 « 45 deg with the motional direction, resulting in a 
Lamb-Dicke parameter r] = k x x o ^ 0.044. Each ex¬ 
perimental sequence begins by initializing the trapped- 
ion oscillator close to its ground state using a com¬ 
bination of Doppler and electromagnetically-induced- 
transparency cooling [25]. We measure a typical mean 
excitation after cooling to be n t h = 0.20(6). Error bars 
throughout this paper are given as standard error of the 
mean. 

In a first set of experiments we characterized our con¬ 
trol using two successive displacements of the potential 
well (Fig. [l]:). At t = 0, we switch the center of the poten¬ 
tial from 0 to Xd, resulting in the ion being excited and 
subsequently oscillating in the displaced potential well. 
At t = At we reverse this change, bringing the potential 
back to its original position. In the ideal case, the ion is 
then in a coherent state of size \a\ given by equation ([!]). 
In order to read out the oscillator state, we then perform 
an analysis using an optical probe pulse of duration t p 
which is resonant with the 5 th sideband. This is followed 
by a detection of the ion’s internal state using resonance 
fluorescence (Supplementary Information). By repeating 
the experiment a large number of times for each value 
of t p , we obtain an estimate of the probability of finding 
the ion spin down as a function of time. This agrees well 
with the functional form 

^t(^p) = 9 Pn t h,M ( n ) [l T e p COS (ft n?n _|_ s t p )J , (3) 

n>0 

where Pn th ,\a\( n ) is the occupation of the n th number 
state for a displaced thermal state with thermal mean 
quantum number n t h and displacement \a\ (this distri¬ 
bution is given in the Supplementary Information). The 
effects of spin and motional decoherence during the probe 
pulse are accounted for by an exponential decay at rate 
T [12]. Experiments were performed for a range of values 
of At, in each case taking data using the first sideband 
s = 1 for 0 < t p < 100 ps. For a single value of At, 
the data do not allow the extraction of |oj, Qq and the 
decay parameter T independently. We therefore probe 
first the non-displaced state directly after cooling. Fit¬ 
ting using equation © with \a\ = 0 and n t h = 0.2, we 
determine D 0 /(27 t) = 181(1) kHz and T = 2.2(3) ms -1 . 
These are then fixed in fits to the data for each value of 
At, from which we extract the corresponding displace¬ 
ment | a\. The resulting values are plotted in Fig. §»■ 
Also shown is a fitted curve using equation ([l]) , which al¬ 
lows us to determine the oscillation frequency of the ion 
cj m /(27r) = 2.3505(6) MHz and amplitude a o = 5.11(1), 
corresponding to x a = 75.0(1) nm (more details of this 
analysis can be found in the Supplementary Informa¬ 
tion). We observe in the results in Fig. that the ampli¬ 
tude of each oscillation of the ion appears to be different. 
We attribute this to drifts in the laser intensity at the 
ion, which result in shifts in the deduced value of \a\ 
(this drift is illustrated by the bounding curves in Fig. 
pk, see Supplementary Information for further details). 
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FIG. 2: Coherent-state generation: a) The states are prepared in the original potential by executing the experimental 
sequence in Fig. [TJ:. The solid orange line is a fit to the data using equation 0- The fluctuations in amplitude can be 
explained by drifts in the laser power over the time taken to acquire the data (shaded area, see Supplementary Information), b) 
Sub-nano-second resolution scan around one complete motional period (A t ~ 27r/cL; m ). This measurement was performed with 
uj m /(27r) ~ 2.53 MHz and corresponds to ao ~ 85. c) Measured probability of finding the spin in the |f) state as a function 
of the probe time using the s — 1 sideband. Filled symbols correspond to the initially-prepared state and open symbols 
after approximately one motional period. This measurement was performed with cj m /( 2tt) « 2.52 MHz and nth ~ 0.15, and 
corresponds to ao ~ 85. Every data point consists of 1000 measurements. Error bars are assigned assuming quantum projection 
noise according to Ref. [24]. 


In order to verify that the switching between poten¬ 
tials does not produce additional effects beyond a simple 
displacement, we examine our ability to catch the ion in 
the ground state after a single cycle of oscillation. In an 
independent measurement we examined data for the re¬ 
gion around At = 2tt /cj m for u; m /(27r) ~ 2.53 MHz and 
an intermediate coherent state of ao ~ 85. To faithfully 
return these states to the origin after a single trap cycle 
we need to trigger the experimental sequence from the 
phase of the radio-frequency trap drive. This can be at¬ 
tributed to a pseudopotential gradient along the axis of 
our trap [26]. The results are presented in Fig. and 
show \a\ =0.05(10) for At « 2ir/u m . Figure [2]: provides 
a direct comparison between data for a non-displaced ion 
with n t h = 0.15(3), and one with the same starting tem¬ 
perature which has been displaced to ao ~ 85 and re¬ 
turned one cycle later. A fit to the data for the latter 
gives fit h = 0.21(3) and consistent values for Qq and T. 

Having characterized our control, we perform a sec¬ 
ond set of experiments to map out the interaction be¬ 
tween the ion and the light field as a function of the 
motional excitation. In contrast to the work described 
above, we send the probe pulse to the ion while the po¬ 
tential is centered at x&. The potential well is displaced 
back to x = 0 after the coherent probe pulse has been 
applied, timing this second displacement to happen once 
the ion has completed an integer number of oscillation 


periods. In this way, we prevent the high oscillator ex¬ 
citation from affecting the internal state detection. Ex¬ 
perimental measurements of P± are shown in Fig. [3^i for 
sidebands between 8 = 0 and s = 5 and displacements 
up to xo ~ 1.5 pan, which corresponds to ao ~ 100 and 
n ~ 10, 000 at the trap frequency of 2.35 MHz. For each 
value of xo we probe P ^ using laser-pulse durations be¬ 
tween 1.4 and 80 qs. To analyze these measurements, 
we start by determining Qo/(2tt) = 204.3(1) kHz and 
n = 0.21(8) from measurements with the non-displaced 
state using the 8 = 0 and 8 = 1 sidebands respectively. 
We then calibrate the displacement sizes using the data 
taken for the carrier transition 8 = 0, fixing T = 0 
and floating only \a\ in equation pi) (for further de¬ 
tails, see Supplementary Information)! These calibrated 
values were then fixed, and data for each of the other 
sidebands was fitted as a full set using a model includ¬ 
ing motional-excitation-dependent AC Stark shifts due 
to nearby atomic transitions and a systematic detuning 
to account for miscalibration (Supplementary Informa¬ 
tion). We find that the AC Stark shifts become relevant 
for 8 = 4 and s = 5 at high motional excitations, be¬ 
cause the frequencies of these transitions are close to the 
seventh and sixth red motional sidebands respectively of 
an atomic transition which is 25.74 MHz away from the 
driven transition. Two-dimensional plots showing the fit 
results are shown in Fig. It below the respective dataset. 
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FIG. 3: Light-atom interaction: a) Measurement of the probability of finding the spin in the |^) state as a function of 
the interaction time with light near-resonant with the s th motional sideband and of the trap-potential displacement x&- Our 
computer-control system imposes t p > 1.4 ps, so for high Rabi frequencies P± can be close to zero even for the first point, b) P± 
obtained using equation with the parameters determined from fits, c) Mean Rabi frequencies for displaced thermal states 
as a function of the potential displacement. The data points are given by the center of the Lorentzian curves to which we fit 
P^(U), the Fourier transform of Pi(t p ). Most error bars are hidden behind the points. The solid lines result from the theory 
model including AC Stark shifts and a detuning (5 0 ff obtained from fitting the data for a given s. 


In order to verify the predicted form of the matrix el¬ 
ements in equation ([ 2 ]), we take a Fourier transform of 
P^(tp) for each value of Xd- This gives us the Rabi fre¬ 
quency distribution, which is peaked due to the limited 
range of number states which contribute significantly to 
the coherent state (Supplementary Information). Fig. 
& shows the mean Rabi frequency obtained from each 
Fourier transform versus the amplitude of the state. We 
overlay this with theoretical curves using equation © 
which include modifications for the offset detunings and 
AC Stark shifts, showing good agreement between ex¬ 
periment and theory over the full range of motional oc¬ 
cupations. The minima in Rabi frequencies observed for 
each set of sideband data are separated by between 515 
and 525 nm. These arise due to the large amplitude of 
the ion’s motion, which spatially samples multiple wave¬ 
lengths of the light. The separation between minima is 
expected to approach the effective wavelength of the light 
as projected on the trap axis for x^ A. In our case this 
value is A/(2cos$) ~ 515 nm (for more details see Sup¬ 
plementary Information). We note that though the 5 = 5 
sideband would take 2 minutes to invert the spin if driven 
in the ground state, the modulation of the light due to 


the large ion oscillation allows us to do this in less than 

10 ps. 

The work presented has a number of possible applica¬ 
tions for quantum information science. It may be used 
to speed-up the transport of trapped-ion qubits in scal¬ 
able quantum information processing [[20ff22] . where the 
transport would proceed by switching the position of the 
potential well to two positions x<\ and 2x d at times which 
are separated by At = 7r/uj m [8]. This would result in the 
ion being “caught” in the ground state of the potential 
centered at 2xa, with the transport taking half a period 
of the ion’s oscillation. We expect that large squeezed 
states could be achieved through sudden trap frequency 
changes, in this case separated by delays of one quarter of 
the trap period [8]. Bang-bang control routines based on 
sudden changes to the trap frequency have been proposed 
for suppressing motional decoherence effects m- Using 
ion chains, the ability to change the trapping potential 
on timescales which are fast compared to the character¬ 
istic frequencies of the interaction between ions has been 
proposed as a method for generating high levels of con¬ 
tinuous variable entanglement by means of the Coulomb 
interaction mmm- 
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FIG. 4: Surface-electrode trap: Linear Paul trap in final 
fabrication steps. Ions are trapped in the position denoted 
with a star, approximately 50 pm from the trap surface. 


I. SUPPLEMENTARY MATERIAL 
A. Experimental setup 

Our experiments are performed using a surface- 
electrode linear radio-frequency (rf) trap in a 5-wire 
asymmetric configuration undo] (Fig. [4]). The pseudo¬ 
potential null line is ~ 50 pm above the trap chip. The 
main trapping zone, denoted with a star in the figure, 
is between the pair of electrodes e5, closer to e5r than 
e51 due to the asymmetry in the width of the rf elec¬ 
trodes. Electrodes e2, e8 and em are each connected to 
the output of a Single-Pole-Triple-Throw (SPTT) switch 
mounted on the Cryo-Electronics Board (CEB, Fig. [5J>) 
and within 3 cm of the trap chip. Voltage offsets of less 
than a volt at these five electrodes suffice to carry out our 
experiments. The trap is driven with an rf amplitude of 
~ 100 V and a frequency of ~ 93 MHz, leading to radial 
secular motion at ^ 4 and « 7 MHz at an axial frequency 
of cj m /(27r) ~ 2.5 MHz. 

The trap is placed in a sealed chamber cooled down 
to ~4K. Cryogenic setups yield improved vacuum over 
room-temperature experiments EQ, resulting in longer 
ion lifetimes. We take advantage of the low outgassing at 
cryogenic temperatures to use standard Printed-Circuit- 
Board (PCB) assemblies for in-vacuum electronics. This 
reduces the technical effort compared to preparing ultra- 
high-vacuum-compatible electronics for room tempera¬ 
ture operation [ [32] . 

The CEB holds the fast switching electronics, 30 low- 
pass filters for the analog input lines, and a track for 
guiding the trap rf drive from a quarter-wave helical res¬ 
onator to the trap chip. For voltage switching we use 
a commercial CMOS integrated circuit (CD74HC4066M, 
from Texas Instruments), which implements four bilat¬ 
eral Single-Pole-Single-Throw (SPST) switches. The ad¬ 
ditional circuitry to implement a SPTT switch is shown 
in Fig. if- The CEB includes five such copies, one per 
switchable electrode. Measured parameters of this circuit 


TABLE I: Summary of the measured performance of the 
switching-electronics at 300 and 4 K. Off and on resistance 
are measured directly at the switch. The power consumption 
was measured at a switching rate of 1 MHz. The crosstalk 
indicated is the fraction of power from the digital side leaking 
to the analog output, and is measured around 1.5 MHz. The 
rise/fall times are given for the 10%-90% transition of a 9 V 
step. 



300 K 

4 K 

Off resistance (MH) 

> 10 

> 10 

On resistance (H) 

12 

3.5 

Power consumption (mW) 

8 

10 

Crosstalk (dB) 

-50 

-50 

Rise/Fall time (ns) 

13 

< 5 


are shown in Tab.[I| Resistors and capacitors are attached 
in parallel pairs in order to have a spare connection in 
the event of solder-point damage during cool-down. To 
ensure cryo-compatibility, all resistors are thin-film and 
capacitors are from the Panasonic ECHU(X) series. Oth¬ 
erwise, the PCB design and soldering techniques do not 
differ from those used to prepare standard non-vacuum 
boards. 

The digital pulses for controlling the switches are pro¬ 
duced by room-temperature electronics based on a mul¬ 
tichannel delay/pulse generator (P400 from Highland 
Technology). Sequences of pulses are triggered by a single 
TTL line locked to the phase of the rf drive and generated 
by a Field-Programmable-Gate-Array (FPGA), which we 
also use for the generation of laser-pulse sequences. 


B. Cooling and Detection 

The experimental sequences are depicted in Fig. (TJ^-d. 
Ground-state cooling is achieved in two stages. The first 
involves Doppler cooling into the Lamb-Dicke regime by 
applying a laser 10 MHz detuned below the Si/ 2 ++ P 1/2 
transition at 397 nm for 500 qs. We then follow this 
with electromagnetically-induced-transparency cooling, 
involving the application of two 397 nm laser fields po¬ 
larized so that they drive simultaneously the a~ and 
7r transitions between the Si/ 2 and Pjy 2 Zeeman sub- 
levels. With suitable combinations of laser powers and 
frequency detunings, this allows us to cool the ion close 
to the ground state of motion [55]. In our setup we reach 
a steady state at n t h ~ 0.2 after 100 qs, limited by the 
heating rate of the ion, which is 1 — 2 quanta/ms from 
the quantum ground state. After cooling, the state is ini¬ 
tialized to |^) by optically pumping with a -polarized 
resonant 397 nm beam. In all these pulses, 866 nm light 
is used to repump the population from D 3 / 2 back to 
Pi/ 2 . Coherent operations between the Si/ 2 and D 5 / 2 
states are performed with a narrow-linewidth 729 nm 
laser. These couple the motional state to the spin sys- 









FIG. 5: Cryogenic electronics: a) Realization of a SPTT switch for a switchable control electrode, r and c form low-pass 
filters for the supply lines of the CMOS switches (V+, VI), with a cutoff frequency of « 7 kHz. R and C form low-pass filters for 
the input DC voltages (14), with a cutoff frequency of « 7 Hz. As usual in rf traps, the capacitors have the additional function 
of rf-grounding the control electrodes. R pro t and R pu n are over-current protection and pull-down resistors for the digital control 
lines of the SPSTs. V a is passed on to the electrode by closing the SPST with a digital high/low level through the control line 
Va- b) Picture of the CEB. The circuitry relative to one of five switches is highlighted in green, a group of seven RC filters in 
orange, and the position of the rf track (inside the PCB) in yellow. The CEB is mounted on a copper structure for thermal 
contact to the cold plate. 


tem, which is subsequently read out by state-dependent 
fluorescence with resonant 397 nm and 866 nm light. In a 
detection time interval of 500 ps, we detect a mean num¬ 
ber of 2 photons for an ion initially prepared in |t), and 
25 photons if the ion is prepared in |j,). 


C. Displaced thermal states 


For a thermal state with average phonon number n t h 
displaced by a , the occupation of the n th number state 
is 33 


Pn th ,\a\(n) = e '“I 2 J2 


n th 


m> 0 


(nth + l ) n+1 


i i | i2 (n—m) 
7 \m\\cv\ v J 


(-1 y la 


2j 


Y~ 

j\(m — j)\(n — m + j)\ 


(4) 


D. Light-atom interaction 


The coupling between the internal electronic state and 
the ion’s motion is described by the Hamiltonian [26] 


frl _ - i ?7 ( ae~ iuJrntp + d t e ia ’ mtp ) — i<5t p 

12 c 2 e e 


+ H.c. 


(5) 


Here, d + = |t)(i| is the spin-flip operator producing 
transitions from ||) to |t) at the Rabi frequency Do/(2tt), 
the Lamb-Dicke parameter r] = k x x o relates the motional 
ground-state wave-packet size xq = yjh/ (2Mcj m ) to the 
projection of the radiation wave vector k on the trap 
axis, M is the ion’s mass, and S is the laser detuning 
from the carrier transition at frequency ojq. For inte¬ 
ger values of 5, laser light detuned by S ~ <sco> m drives 
near-resonant transitions on the 5 th motional sideband 
|n) |t) |n+s). The corresponding matrix element is 

proportional to \(n + s\e lkxX \n) | and its associated Rabi 
frequency can be evaluated analytically [23] : 


a 


n,n+s 


— ^0 


(n + s|e i?7(at+a) |n> 

= (r? 2 )| , (6) 


with I = max {0, m — n}. 


where, n < (n>) stands for the lesser (greater) of n + s 
and n, and L%(x) are generalized Laguerre polynomials. 
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number state n 

FIG. 6: AC Stark shifts: Energy level corrections due to 
off-resonant coupling to motional sidebands other than the 
one driven, as well as to the carrier and sidebands of the 
\i) ++ \L = 2, J — 5/2, Mj = 3/2) transition at our 3.83 G 
magnetic field when 5 = 0 at n = 0. 


E. Dependence of deduced oscillation amplitude on 
carrier Rabi frequencies 


lustrate this we repeat the fits fixing the Rabi frequency 
to the limiting values Do =•= £Do and floating olq and cj m . 
This yields ao = 5.00(1) for D 0 ± £D 0 , and = 5.22(1) 
for Dq — SSl 0 , with u; m /(27r) = 2.3505(6) MHz in both 
cases. The boundaries of the shaded area in Fig. [2^i are 
obtained by inserting these numbers into equation ([!]). 


F. Motional-state dependent AC Stark shifts 


The model used for calculating the expected in our 
experiment (Fig. [3| includes a detuning from the driven 
motional sideband as well as an AC Stark shift due to 
off-resonant coupling to transitions other than the one 
probed. The latter is given by 


S A c(n,s) 


E 


Q 2 

^ L n,n+s 

2w m (s' - s ) 


( si n>n+e /w_i ^) 2 

~f 4(A_i_).| + w m (s' - s)) 


(7) 


For the first set of experiments we carried out the se¬ 
quence in Fig. [ 1 J 3 and obtained the results in Fig. [2] Here 
we focus on one of the consequences of the slow drifts in 
laser power, which are not included in the model used for 
theoretical calculations. 

From a first fit of a non-displaced state we obtain 
D 0 /(27r) = 181(1) kHz and T = 2.2(3) ms -1 , using equa¬ 
tion © with | a | = 0 and n t h = 0.2. From subsequent 
fits for each value of At we obtain the coherent-state 
sizes plotted in Fig. [2k. Fitting these to equation ([!]), we 
determine = 5.11(1) and co’ m /(27r) = 2.3505(6) MHz. 
However, we note that over the time required to take this 
data the Rabi frequency varies over a range $Do/(27r) ~ 
±3.2 kHz. The value of \a\ determined from fits is 
strongly correlated to the Rabi frequency used, which 
accounts for the fluctuations observed in Fig. nt- To il¬ 


and plotted in Fig. [6] for our experimental parameters. 
The first term arises from coupling to sidebands of the 
spin transition, while the second is due to the fact that 
the 729 nm laser drives off-resonantly the secondary tran¬ 
sition ||) \L = 2, J = 5/2, Mj = 3/2). The coupling 
strength to the secondary transition is reduced relative 
to that of the resonantly driven transition by a factor 
W_ i_^| = y/b. At our magnetic field of ~ 3.83 G, 
the frequency gap between both carrier transitions is 
A_i^|/(27r) « 25.74 MHz. This is close to an inte¬ 
ger times the motional frequency, causing the s = 4 (5) 
sideband of the main transition to be almost resonant 
with the s = —7 (—6) sideband of the secondary tran¬ 
sition, and therefore strongly AC Stark shifting the ||) 
state. For a total detuning 5 t ot(^? s ) = ^off(^) ±£ac(?T s), 
equation © can be re-written as 


P^th, M = 1 - E^,|a|( n )^2- 


^n,n+s 


/ r>2 



( 8 ) 


where we have ignored the exponential-decay term due 
to decoherence, since our probe times go up to 80 ps and 
we found previously r^2ms _1 . 

G. Light-matter-interaction data analysis 

For this experiment we carried out the sequence in Fig. 
Eh The results are shown in Fig. It- To analyze the 


data and compare it to the theoretical model, we first 
calibrated the displacements with the data obtained for 
the carrier (5 = 0 ). We fitted the measurement results for 
each value of x d to equation © while fixing Do and n to 
values determined previously by driving Rabi oscillations 
on the 5 = 0 and 5 = 1 transitions, respectively. This 
yielded a coherent-state size which we converted into dis¬ 
placements according to ao = Xd/(2xo). For the rest of 
the sidebands we used the calibrated displacements and 
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TABLE II: Fit results for the light-matter-interaction exper¬ 
iment. The probe-pulse frequencies for sidebands s = 0 to 2 
were pre-calibrated with non-displaced states. Higher-order 
sidebands are too weak for a direct measurement with non- 
displaced states, so probe-pulse frequencies for s > 2 where 
estimated from the lower sideband frequencies. This explains 
the large systematic shifts in 5 0 s resulting from the fits. 


s 

Oo/(27r) (kHz) 

<5off/( 2tt) (kHz) 

0 

205 

0 

1 

211 

4 

2 

218 

-4 

3 

224 

-12 

4 

226 

-23 

5 

226 

-32 



Rabi frequency (kHz) 


FIG. 7: Rabi-frequency spectrum: The data points and 
curve are obtained as explained in the text. This plot is for 
the carrier transition (s = 0) and a displacement xa ~ 660 nm 
(a 0 « 44). 


fitted the measurement data to find and £ 0 ff for each 
individual sideband. The best fit parameters are given 
in Tab. [II] and yield the results in Fig. Ilf- c. 

To analyze the evolution of the mean Rabi frequency 
for a given driven sideband s and as a function of the state 
size (Fig. [3^), we calculate the discrete Fourier transform 
of the measured P|(t p ) and propagate the shot-noise un¬ 
certainties to the frequency domain according to [ [ 34 ] . 
The number-state distribution Pn th ,\a\( n ) for low values 
of n t h is narrow compared to the features which result 
from equation so the frequency spectrum shows a 
single peak at the mean Rabi frequency (Fig. [7]). We fit 
these to symmetric Lorentzian functions, from which we 
determine the center frequencies given in the plots. We 
compare this to the theoretical mean Rabi frequency for 
a displaced thermal state, given by 

^th = y^Pn th ,|q| ( n )\J ^n,n+s ^tof (9) 

n 

The separation between minima in mean Rabi fre¬ 
quency depends on the radiation wavelength A. For large 
values of n, the Laguerre polynomials in equation |6| can 
be approximated by Bessel functions J a as [ 35 ] 



The zeros in the right-hand side correspond to zeros 
of the Bessel functions, whose separation tends to i r 
for large arguments. This implies that the separa¬ 
tion between minima in mean Rabi frequencies tends to 
A/(2 cos 0) for large values of x&. 








